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1 Introduction 

The strict positive realness (SPR) of transfer functions is an important performance specification, 
and plays a critical role in various helds such as absolute stability/hyperstability theory [15, 20], 
passivity analysis [11], quadratic optimal control [3] and adaptive system theory [16]. In recent years, 
stimulated by the parametrization method in robust stability analysis [1, 4, 6], the study of robust 
strictly positive real systems has received much attention, and great progress has beem made [2, 5] [7]- 
[10] [12]-[14] [17]-[19] [21]-[35]. However, most results belong to the category of robust SPR analysis. 
Valuable results in robust SPR synthesis are rare. The following fundamental problem is still open 
[14, 18, 22, 23] [25]-[28] [30] [31]-[34]: 

Suppose a{s) and b(s) are two n-th order Hurwitz polynomials, does there exist, and how to find a 
(fixed) polynomial c(s) such that c(s)/a(s) and c{s)/b{s) are both SPR? 

By the definition of SPR, it is easy to know that the Hurwitz stability of the convex combination 
of a{s) and 6(s) is necessary for the existence of polynomial c(s) such that c{s)/a(s) and c(s)/6(s) are 
both SPR. In [13, 14, 19], it was proved that, if a(s) and b(s) have the same even (or odd) parts, such a 
polynomial c(s) always exists; In [2, 13, 14, 17, 22, 23, 25, 27, 30, 31], it was proved that, if n < 4 and 
a(s),6(s) S K (K is a stable interval polynomial set), such a polynomial c(s) always exists; Recent 
results show that [22, 23, 25, 27, 28, 30] [32[-[34], if n < 6 and a(s) and b(s) are the two endpoints of 
the convex combination of stable polynomials, such a polynomial c(s) always exists. Some sufficient 
condition for robust SPR synthesis are presented in [2, 5, 10, 17, 22, 23, 27], especially, the design 
method in [22, 23] is numerically efficient for high-order polynomial segments or interval polynomials, 
and the derived conditions are necessary and sufficient for low-order polynomial segments or interval 
polynomials. 

This paper presents a constructive proof of the following statement: for any two n-th order poly¬ 
nomials a(s) and b{s), the Hurwitz stability of their convex combination is necessary and sufficient for 
the existence of a polynomial c(s) such that c(s)/a(s) and c(s)/b(s) are both SPR. This also shows 
that the conditions given in [22, 23] are also necessary and sufficient and the above open problem 
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admits a positive answer. Some previously obtained SPR synthesis results for low-order polynomial 
segments [22, 23, 25, 27, 28, 30] [32]-[34] become special cases of our main result in this paper. 


2 Main Results 

In this paper, P" stands for the set of n-th order polynomials of s with real coefficients, R stands for 
the field of real numbers, d{p) stands for the order of polynomial p{-), and iJ" C P" stands for the 
set of n-th order Hurwitz stable polynomials with real coefficients. 

In the following definition, p{-) G P™, q{-) G P", f{s) = p{s)/q{s) is a rational function. 

Definition 1 [29] /(s) is said to be strictly positive real(SPR), if 

(i) d{p) = d{q); 

(ii) /(s) is analytic in Re[s] > 0, (namely, q{-) G P" ); 

(hi) Re[/(jw)] >0, Vw G R. 

If /(s) = p{s)/q{s) is proper, it is easy to get the following property: 

Property 1 [9] If /(s) = p{s)/q{s) is a proper rational function, q{s) G P", and Vw G 
P, Re[/(jw)] > 0, then p{s) G P” U P”-b 

The following theorem is the main result of this paper: 

Theorem 1 Suppose a{s) = s” -|- + • ■ • + Un G P", b{s) = s” -I- + •■• + &« S P", 

the necessary and sufficient condition for the existence of a polynomial c(s) such that c(s)/a(s) and 
c(s)/6(s) are both Strict Positive Real, is 

A6(s) + {1- A)a(s) G P”, A G [0,1]. 

The statement is obviously true for the cases when n = I or n = 2. We will prove it for the case 
when n > 3. 

Since SPR transfer functions enjoy convexity property, by Property I, we get the necessary part 
of the theorem. 

To prove sufficiency, we must first introduce some lemmas. 

Lemma 1 Suppose a{s) = s” -I- + • ■ • + Un G P", then, for every k G {I, 2, • • •, n — 2}, 

the following quadratic curve is an ellipse in the first quadrant (i.e., Xi > 0,i = I, 2, • • •, n — I) of the 
space {xi,X 2 , • • •, Xn-i) 

{ Cfc+i - 4CfcCfc+2 = 0, 

Cl = 0, 

I G {1, 2, • * •, Tz}, I ^ kjk -\- \^k -\- ‘1^ 

^When n = 3 , the ellipse equation is (see [ 27 , 28 ] for details): 

(a2Xi — aiX2 - as)^ - 4 (ai — xi)(a 3 a: 2 ) = 0 . 

When n = 4 , the two ellipse equations are (see [ 25 ] [ 30 ]-[ 32 ] for details): 

{ (<12X1 -t- X3 — aiX2 — as)^ — 4 (ai — xi)(a3a;2 — <12x3 — 04x1) = 0 , 

04x3 = 0, 

{ (a 3 X 2 — a2X3 — a 4 Xi)^ — 4 (a 2 Xi -|- X3 — aia;2 — a 3 )a 4 ,X 3 = 0, 
ai — xi = 0. 

When n = 5 , the three ellipse equations are (see [ 33 ] for details): 

f (a 2 Xi -j- X 3 — aiX 2 — a3)^ — 4 (ai — a;i)(a5 + a 3 X 2 -t aia;4 — 02X3 — 04X4) = 0 , 

\ a 4 X 3 — 03x4 — a 5 X 2 = 0, a 5 X 4 = 0, 

/ (as + a 3 X 2 + aiX 4 — a2X3 — 04X4)^ — 4(02X4 + X3 — 04X2 — a3)(a4X3 — 03X4 — 05X2) = 0 , 

\ 04 — X4 = 0 , 05x4 = 0 , 

f (04x3 — 03x4 — 05x2)^ — 4(05 + 03x2 -I- 04x4 — 02x3 — 04x4)05x4 = 0 , 

\ 04 — X4 = 0 , 02 X 4 -I- X3 — 04x2 — 03 = 0 . 
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and this ellipse is tangent with the line 


f Q = 0, 


and the line 


f Q = 0, 

\ I € {1,2,-,n},l ^ k,k + I, 


respectively, where ci := ^ ajX 2 i-j-i, ^ = 1, 2, • • •, n, oq = 1, xq = 1, Oi = 0 if i < 0 or i > n, 

j=o 

and Xi = 0 if i<0 or i>n—1. 

Proof Since a{s) is Hurwitz stable, by using mathematical induction. Lemma 1 is proved by a 
direct calculation. 

For notational simplicity, for a{s) = s” + ais”“^ + - • ■ + a„ G iL”, b{s) = s” + 6is”“^ + - • • + &„ G i?", 
Vfc G {1, 2, • • •, n — 2}, denote 

4+i - 4cfcCfc+2 < 0, 

Cl = 0,/ G + 

and 

^Ik ■= {ixi,X 2 , - ■ ■ ,Xn-l)\ dfc+i - 4:dkdk+2 < 0, 

di = 0,1 G {l,2,---,n},/ ^ k,k + l,fc + 2}, 

n n 

where c; := X] d; := X) ^ = l,2,---,n, Oq = l,6o = l,a^o = 1, 

j=0 j=0 

tti = 0 and = 0 if i < 0 or t > n, and a;i=0ift<0ori>n — 1. 

In what follows, {A, B) stands for the set of points in the line segment connecting the point A and 
the point B in the space {xi,X 2 , • ■ •, Xn-i), not including the endpoints A and B. Denote 

'■= {(xi,X 2 ,---,X„_i)l (xi,X 2 ,---,X„_i) G 

VA, GL!“„zG{l,2,--.,n-2}} 

and 

:= {(xi, X2, * * * 5 ^n—1)| (^1) ^2; ' ‘ ' 7 ^n — l) ^ Ui^l z<j<n— 

vs, GD^,iG{l,2,--.,n-2"}}. 


Lemma 2 Suppose a(s) = s” + ois” ^ + ■ ■ ■ + a„ G id", b(s) = s” + &is” ^ + ■ ■ ■ + b„ G id", 

if 0“ n 7 ^ (j), take {xi,X2, • • • , a;„_i) G fi" n and let c(s) := s"“^ + (xi — £)s”“^ + X 2 s"“^ + 

c(s) c(s) 

—^ and 77 -^, we nave 
a(a) b{s) 

Vw G R, Re[4^] > 0 and Re[^^] > 0. 

a(jw) b{juj) 

Proof Suppose {xi,X 2 , • • •, a;„_i) G D let c(s) := s”“^ + (xi — £:)s"“^ + X 2 s"‘~^ + • • • + 
Xn- 2 S + {xn-1 + ff), £ > 0, £ Sufficiently small. 

When n = 6 , the four ellipse equations are (see [ 34 ] for details): 

r (a2Xi + X3 — aiX2 - as)^ — 4 (ai - 3 ;i)(a 5 + a 3 X 2 + aia;4 — X5 — 02X3 - a4Xi) = 0 , 

\ aexi + 04x3 + 02X5 — 03x4 — 05x2 = 0, 05x4 — 04x5 — 05x3 = 0, 05x5 = 0, 

(05 + 03X2 + 04X4 — X5 — 02X3 — 04X1)^ — 4(02X1 + X3 — 01X2 — 03)(06Xi + 04x3 + 02X5 — 03X4 — 05x2) = 0 , 
oi — XI = 0, 05x4 — 04x5 — 05x3 = 0, 05x5 = 0, 

(ogxi + 04x3 + 02x5 — 03x4 — 05x2)^ — 4(05 + 03x2 + 04x4 — X5 — 02x3 — a4Xi)(a5X4 — 04X5 — 05X3) = 0 
01 — XI = 0, 02X1 + X3 — 01x2 — 03 = 0, 05x5 = 0, 

(05x4 - 04x5 - 05x3)2 - 4(05x1 + 04x3 + 02x5 - 03x4 - 05x2)05x5 = 0, 

01 — XI = 0, 02x1 + X3 — 01x2 — 03 = 0, 05 + 03x2 + 01x4 — X5 — 02x3 — 04x1 = 0. 


• • • + Xn-2S + (xn-1 + s) {s is a Sufficiently small positive number), then for 
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Vw G R, consider 


Re[ 


c(j^) j 

a(jw) 


a{juj) 


r[ciw^(” ^) + • • • + c„_ia;^ + c„] 


H-Re[ 




+ 


aO'u;) 

H-h Cn-lUj'^ + Cn] 




where q := ^ (—l)^+'’ajX 2 i-j-i, ? = 1, 2, • • •, n, ao = 1) a^o = 1) = 0 if i < 0 or i > n, and Xj = 0 if 

i=o 

i<0 ori>n — 1, and c(oj^) is a real polynomial with order not greater than 2 (n — 2 ). 

In order to prove that Vw G R, Re[ '"^/^| 1 > 0, let t = to'^, we only need to prove that, for any 

a{juj) 

£ > 0,e sufficiently small, the following polynomial fi(t) satisfies 

fl{t) := + C 2 t”“^ H-h Cn-lt + c„ 

- c{t)) >0, Vt G [0, +CX)). 


Since (xi, X 2 , • • •, x„_i) G by the definition of 12“, it is easy to know that 

gi{t) := + C 2 t"“^ H-h c„_it + c„ > 0, Vt G (0, +oo). 

Moreover, we obviously have /i(0) > 0, and for any £ > 0, when t is a sufficiently large or sufficiently 
small positive number, we have /i(t) > 0 , namely, there exist 0 < ti < ^2 such that, for all e > 0 , 
t G [0,ti] U [^ 2 , +oo), we have /i(t) > 0. 

Denote 

Ml = inf 5 i(t), 
tG[tl,t2] 

Ni= sup |r-^-c(t))|. 
tG[ti,t2] 

Ml 

Then Mi > 0 and Ni > 0. Choosing 0 < £ < , by a direct calculation, we have 

Ni 

fi{t) := + C 2 t”“^ H-|-c„_it + c„ 

+£(t”“^ - c(t)) >0, Vt G [0, +CX)). 

Namely, 

Vw G i?,Re[4^] > 0. 

Similarly, since (xi, X 2 , • • •, x„_i) G there exist Q < such that, for all £ > 0, t G 

[ 0 ,^ 3 ] U [^ 4 , +00), we have / 2 (f) > 0,where 

72(0 dit" ^ + ^ 2 ^" ^ + • ■ • + dn-it + dn 
+£(t”“^ - d{t)), 


di :='^{-iy+^bjX 2 i-j-i, I = 1,2, - ■ ■ ,n, 

3=0 

where 60 = !> a^o = 1 ) = 0 if i < 0 or i > n, and x^ = 0 if i < 0 or i > n — 1 , and d{u)‘^) is a real 

polynomial with order not greater than 2 (n — 2 ) which is determined by the following equation: 


Re[ 


-e{jujy 


+ £1 




(-g) 

Hj^)V 




dV))- 
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Denote 

92(t) := dit" ^ + fi2t” ^ + • ■ • + dn-it + dn, 
M 2 = inf g 2 {t), 

^€[^ 3 ,^ 4 ] 

N 2 = sup |t”“^ - (i(t))|. 

te[i3.t4] 

M 2 

Then M 2 > 0 and N 2 > 0. Choosing 0 < e < -, we have 

N 2 


Voj G R, Re[ 




> 0 . 


Ml M 2 

Thus, by choosing 0 < e < min{——, ——}, Lemma 2 is proved. 

A^i N 2 


Lemma 3 Suppose a{s) = s” + ms” ^ + • ■ • + a„ G iL", b{s) = s” + 6is” ^ + • • ■ + G iL", if 
\h{s) + (1 — A)a(s) G iL", A G [0,1], then n 12^ ^ 

Proof If VA G [0,1], aA(s) := \b{s) + (1 — A)a(s) G iL", by Lemma 1, for any A G [0,1], 
■■■ ,n — 2, are n — 2 ellipses in the first quadrant of the i?" ^ space (xi, a; 2 , • • •, Xn-i)- 


VA G [0,1], denote 


^ ^ ■ — { (xi , 2^2) ' ■ ' ) 2;^— 1 ) I ( 2^1 ,2^2, ■ ' ■ , Xfi — l) G Ui=iy<j<n—2 (^Ai, AIaj ) , 

VAa* G G {l,2,---,n-”2}} 


Apparently, when A changes continuously from 0 to 1, will change continuously from to 
11 ^, and will change continuously from 11 “^, to Dg^., k = l, 2 ,---,n — 2 . 

Now assume 11“ n = </>, by the definitions of 11“ and and Lemma 1, 3ui > 0, U 2 > 0, ui 7 ^ 
ai,Mi 7 ^ 61 , and 3fc G {l,2,---,n — 2}, such that the following hyperplane L in the space 

(xi, X2 , ' * ‘ , 2 ^n—1) 


L : 


Xl 


X2 

— + ••• + 


2^n—1 


= 1 


Ui U2 Un-1 


separates 11“ and D^, meanwhile, L is tangent with 11“]^, 11“2, • • •, ^e{n- 2 ) ^efc simultaneously (or 
tangent with llg^, llg 2 , • • •, ^^e(„_ 2 ) and 11“^ simultaneously). 

Without loss of generality, suppose that L is tangent with 11“;^, 11“2, • • •, ^e(n- 2 ) and simulta¬ 
neously. 

In what follows, the notation [x] stands for the largest integer that is smaller than or equal to the 
real number x, and {y)^ stands for the remainder of the nonnegative integer y divided by the positive 
integer z 

Since L is tangent with 11“^, 11“2, • • •, ^e(„_ 2 ) and DC simultaneously, note that a(s) is Hurwitz sta¬ 
ble and Ml > 0 , Ml 7 ^ ui, M 2 > 0 , using mathematical induction, by a lengthy calculation, we know that 
the necessary and sufficient condition for L being tangent with D“;^, D“ 2 , • • •, f^e(„_ 2 ) simultaneously 

®For example, [1.5] = 1, [0.5] = 0, [—1.5] = —2, and ( 0)2 = 0, ( 1)2 = 1, ( 11)3 = 2. 
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IS 


and 


? —1“ 1 r ^ 1 r ^ 1 

” [^] (i+n [ 9 ] - [ 9 ] 

^(-1) 2 ^ =Q 

2=0 


\Lil] . [-] 

% = (-!) 2 , j = 3,4, •••,n-1, 


( 1 ) 


— - u-]^ u -2 , 9 — u,-i, •■•,/(,— X, (2) 

where ao = 1. 

[^] O') 

Since Uj = (—1) 2 ^^ 2 ^ > J = 3,4, • • •, n— 1, L is tangent with , by a direct calculation, 

we have 

[^] o+iN [ 9 ] - [ 9 ] 

^(-1) 2 2^0 (3) 

i=0 

where bo = 1. 

From (1) and (3), we obviously have VA S [0,1], 


r* + 1, 


^(-1) 2 OAiUi 


2=0 


( 4 ) 


i=0 


where uai := Oj + X{bi — a^), uq = 1,5o = 1, t = 0,1, 2, • • •, n. (4) and (2) show that L is also tangent 
with (VA S [0,1]), but L separates n“r and fl^r, and when A changes continuously from 0 to 1, 
will change continuously from f2“~ to Xl^r, which is obviously impossible. This completes the proof. 

6rC GK 

Lemma 4 Suppose a(s) = s" + Ois" ^ + • • • + ttn S 7J”, 6(s) — s" + 6is" ^ + ■ ■ ■ + V i?", c(s) — 


c(jw) 


rC(jw), 


s” ^ + xis” ^ +-h Xn- 1 , if Vw S R, Re[——-] > 0 and Re [--—-] > 0, take 

a{juj) b{juj) 

c (s) := c(s) + 6 ■ h{s), S > 0, 6 sufficiently small. 


c ( 5 ) c ( 5 ) 

(where h{s) is an arbitrarily given monic n-th order polynomial), then —and -rr-^ are both 

a(s) b{s) 

strictly positive real. 

Proof Obviously, 9(c) = 9(a) = n, namely, c(s) and a(s) have the same order n. Since a(s) S i4", 

there exists wi > 0 such that, for all \ co |> wi, we have Re(^ 7 ^^) > 0. 
_ a(ju;) 

^"When n = 3 , we have (see [ 27 , 28 ] for details): 

U 1 U 2 — CL 1 U 2 — 0)2^1 + as = 0. 

When n = 4 , we have (see [ 25 ] [ 30 ]-[ 32 ] for details): 

— aiii2 — a2UiU2 + a3'U2 + a4iii = 0, 
uz = -UiU 2 . 

When n = 5 , we have (see [ 33 ] for details): 

— aia2 — a2UiU2 + C13U2 + a4'Ui — as = 0, 

-0.3 = —ai'U2, " 0.4 = — 

When n = 6, we have (see [ 34 ] for details): 

— aiii^ — a2UiU2 + CLZU 2 + a4'Uia2 — cl^U2 — clqui = 0, 

-0.3 = —'Ui'U2, " 0.4 = —U21 Uz = 'Ui'U^. 
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Denote 


M 3 = inf Re(^^l^), N 3 = sup | Re(^^y^) |. 
\UJ\<UJ1 a{juj) \u\<LJi a{juj) 

M3 

Then M 3 > 0 and > 0. Choosing 0 < £ < , it can be directly verified that 

-/V3 


Re(4^) > 0,Vw e R. 
a{juj) 


0 . 


Similarly, 9(c) = d(b) = n, and there exists W 2 > 0 such that, for all I w |> 102 , we have ) > 

o{3^) 

Denote 


M 4 = inf Re(||^), IV 4 = sup | Re(^|l^) |. 
H<u ;2 b(juj) 

Then M 4 > 0 and iV 4 > 0. Choosing 0 < e < it can be directly verified that 

A'4 


Re( 


c(i^) 

^(jw) 


) 


> 0 , Vw e R. 


Thus, by choosing 0 < e < min{—-^}, Lemma 4 is proved. 

N4 

The sufficiency of Theorem 1 is now proved by simply combining Lemmas 1-4. 


Remark 1 From the proof of Theorem 1, we can see that this paper not only proves the existence, 
but also provides a design method. In fact, based on the main idea of this paper, we have developed 
a geometric algorithm with order reduction for robust SPR synthesis which is very efficient for high 
order polynomial segments [26]. 

Remark 2 The method provided in this paper is constructive, and is insightful and helpful in 
solving more general robust SPR synthesis problems for polynomial polytopes, multilinear families, 
etc.. 


Remark 3 Our main results in this paper can also be extended to discrete-time case. In fact, by 
applying the bilinear transformation, we can transform the unit circle into the left half plane. Hence, 
Theorem 1 can be generalized to discrete-time case. Moreover, in discrete-time case, the order of the 
polynomial obtained by our method is bounded by the order of the given polynomial segment [27, 29]. 


Remark 4 
is also SPR. 


If and are both SPR, it is easy to know that VA G [0,1], -—— ,,, , , 
a(s) bis) > J Aa(s)-f (I - A) 6 (s) 


Remark 5 The stability of polynomial segment can be checked by many efficient methods, e.g., 
eigenvalue method, root locus method, value set method, etc. [1, 4, 6 ]. 


3 Conclusions 

We have constructively proved that, for any two n-th order polynomials a(s) and 6 (s), the Hurwitz 
stability of their convex combination is necessary and sufficient for the existence of a polynomial c(s) 
such that c(s)/a(s) and c(s)/ 6 (s) are both strictly positive real. By using similar method, we can 
also constructively prove the existence of SPR synthesis for low order (n < 4) interval polynomials. 
Namely, when n < 4, the Hurwitz stability of the four Kharitonov vertex polynomials is necessary 
and sufficient for the existence of a fixed polynomial such that the ratio of this polynomial to any 
polynomial in the interval polynomial set is SPR invariant [22, 23, 25, 27, 30, 31]. The SPR synthesis 
problem for high order interval polynomials is currently under investigation. 
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